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Abstract 

Stabiiity, bifurcation properties, and the spatiotemporai behavior of different nonlinear combi- 
nation structures of spiral vortices in the counter rotating Taylor-Couette system are investigated 
by full numerical simulations and by coupled amplitude equation approximations. Stable cross- 
spiral structures with continuously varying content of left and right spiral modes are found. They 
provide a stability transferring connection between the initially stable, axially counter propagating 
wave states of pure spirals and the axially standing waves of so-called ribbons that become stable 
slightly further away from onset of vortex flow. 

PACS numbers: 47.20.-k, 47.54.+r, 47.32.-y, 47.10. +g 
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The combination of two counter propagating plane linear traveling waves (TWs) of equal 
amplitude creates a standing wave (SW). This property is also the reason for the possibility 
of SW structures in a variety of nonlinear pattern forming systems that show an oscillatory 
instability and inversion symmetry in one or more spatial directions Then two symmetry 
degenerate counter propagating TW structures as well as a nonlinear SW solution appear at 
the common oscillatory bifurcation threshold. However, in many systems like, e.g., binary 
mixture convection TW solutions bifurcate subcritically being initially amplitude unstable 
and then become strongly nonlinear on the stable upper solution branch. Moreover, the SW 
solution is often amplitude unstable and/or unstable against phase propagation and TW 
formation. 

Thus, we consider here nonlinear combination states of left-spiral vortex (L-SPI) and 
right-spiral vortex (R-SPI) structures that — being mirror images of each other — travel 
either upwards or downwards in the annulus between two counter rotating cylinders of a 
Taylor- Couette setup . This system offers an easy experimental and numerical access 
(i) to forwards bifurcating initially stable TWs (i.e., L-SPI and R-SPI), iii) to forwards 
bifurcating SWs [called ribbons (RIBs) in the Taylor- Couette literature |2j, |3J]] that become 
stable slightly above threshold, and in particular also (iii) to stable cross-spirals (CR-SPI), 
i.e., combinations of L-SPI and R-SPI with different amplitudes which provide a stability 
transferring connection between TW and SW solution branches. This bifurcation scenario 
does not seem to have been reported so far 

Here, we describe stability, bifurcation properties, and the spatiotemporal behavior of 
the three oscillatory structures (i)-(iii) by two amplitude equations with a coupling that 
contains necessarily also a quintic order term to enable CR-SPI solutions. Their results are 
compared with numerical solutions of the full Navier-Stokes equations that we obtained for 
a system with radius ratio r] = 1/2 by methods described in |9]. 

The vortex structures (i)-(iii) are axially and azimuthally periodic with wave numbers 
k = 2tt/\ and M, respectively, with A = 1.3 and M = 2 throughout this paper. They 
rotate with characteristic constant angular velocities into the same direction as the inner 
cylinder :QJ and thereby they are forced to propagate axially except for RIB vortices which 
only rotate but do not propagate. 

In order to characterize the spatiotemporal properties of the vortices we used also Fourier 
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decompositions in azimuthal and axial direction 

f(r,<p,z,t) = Y,f™Ar,t)e^ +nkz K (1) 

m,n 

For example, L-SPI (R-SPI) flow contains only modes with m = nM (m = —nM) since it 
does not depend on ip, z, t separately but only on the phase combination <pA = Mip + kz — UAt 
{<pB = Mip — kz — ujBt). Here we use A and B to identify L-SPI and R-SPI properties, 
respectively. For SPI that are axial mirror images of each other the modes in oscillate 
according to fm,n( r ,t) = $m,±nM e ^ f±n{ r ) with a common frequency uja = ojb = 

ui(k,M). Thus, the color-coded plots in Fig. GJa,b) of the SPI in the ip — z plane of an 
'unrolled' cylindrical surface are patterns of straight stripes. They both move in positive 
^-direction, i.e., into the direction of rotation of the inner cylinder Q with angular phase 
velocity <pspi = oj/M. The L-SPI propagates upwards with axial phase velocity wA = u/k 
and and its mirror image the R-SPI propagates downwards, Wp h = —u/k . 

On the other hand, a CR-SPI is a kind of combination of a L-SPI and a R-SPI with 
different flow amplitudes and different frequencies. Thus, the CR-SPI spectrum contains 
two frequencies, say, us and Ud, that are in general incommensurate. Within the amplitude 
equation approximation presented below us = {oja + u)b)/2 and ujd = (<^>A ~ ^b)/2. For the 
CR-SPI in Fig. [IJc) the amplitude of the L-SPI content is larger than the R-SPI content 
while uja < ujb- Then the periodic deformation of the L-SPI-dominated stripe pattern in 
Fig. njc) that is caused by the minority R-admixture propagates axially downwards with 
velocity uJo/k while rotating with angular velocity us/M . The CR-SPI in Fig. [D^d) is the 
mirror image of Fig. ^c) with reverse amplitude and frequency relations. For all CR-SPI 
investigated here us ^> \ujd\- RIB states are degenerate CR-SPI states with equal content 
of R-SPI and L-SPI contributions and ur> = which do not propagate axially. 

Fig. Elshows for two different R\ bifurcation diagrams of SPI, CR-SPI, and RIB solutions 
versus R 2 and its reduced distance \i = (R 2 — R^)/\R^\ from the onset, R®, of SPI and 
RIB flow. Here R\ and R 2 are the Reynolds numbers defined by the rotational velocities 
of the inner and outer cylinder, respectively. Order parameters are the squared amplitudes 
(|A| 2 , \B\ 2 ) and frequencies (u)a,wb) of the dominant modes 

u 2>1 (t) = \A\e~ iuJAt , U2,-i(t) = \B\e~ iuJBt (2) 

in the decomposition ((TJ) of the radial velocity u at mid-gap. For all vortex structures 
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investigated here the moduli and frequencies of u 2 ,±i in Eq. (J2J) are constant. 

The L-SPI (A ^ 0, B = 0) shown in Fig.|2]by circles and the RIB solution (A = B) marked 
by diamonds bifurcate at /j = with common linear frequency uj° out of the unstructured 
CCF. However, with increasing \i their frequencies vary differently. The squared amplitudes 
of these two states grow basically linearly with /x, albeit with different slopes. 

Initially, the SPI is stable and the RIB is unstable. But then there appears at larger 
H a stable CR-SPI solution (triangles) which transfers stability from the SPI to the RIB. 
The CR-SPI {A ^ B) bifurcates with B = and finite ujb out of the L-SPI as shown in 
Fig. El (bifurcation out of the symmetry degenerate R-SPI is analogous — just exchange A 
and B). Then \B\ grows and \A\ decreases while u>a increases and ujb becomes smaller until 
the CR-SPI solution branches end with A = B and uja = in the RIB state. The latter 
loses stability outside the plot range of Fig. [21 to amplitude-modulated CR-SPI that are not 
discussed here. 

Spatiotemporal properties, bifurcation, and stability behavior of the SPI, RIB, and CR- 
SPI states close to onset can reasonably well be explained and described within an amplitude- 
equation approach. Therein the fields are represented by the superposition — here written 
down, e.g., for the radial velocity u(r, <p, z, t) — 

u = [A(t)u A (r)e ikz + B(t)u B (r) e - lkz ] e ^~^ + c.c. (3) 

of just the linear critical modes / u^ v B(r)e ± * fc2 e^ M¥ '~ a;0 ^ of L- and R-SPI with slowly varying 
amplitudes A and B, respectively. Here we normalize \ua\ and \ub\ to 1 at mid-gap. There 
the representation (jSJ) is quite sufficient. But closer to the inner cylinder the vortex fields 
contain further axial and azimuthal modes that are generated via nonlinear interactions 
of the critical modes. However, it seems that close to onset they do not influence deci- 
sively the bifurcation and stability behavior that is governed by the modes retained in the 
approximation ©. 

In order to reproduce the bifurcation and stability behavior of the aforementioned vortex 
states including the CR-SPI one needs coupled equations for A and B of at least quintic 
order [l^] • Higher-order terms that are suggested in [^J are not necessary to ensure existence 
of CR-SPI solutions. Symmetry arguments |3] and simplicity considerations discussed below 
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restrict the form of the equations to 

A=[afi+ b\A\ 2 + c\B\ 2 + e (\A\ 2 - \B\ 2 ) \B\ 2 ] A (4a) 
B=[afi + b\B\ 2 + c\A\ 2 + e (\B\ 2 - \A\ 2 ) \A\ 2 ] B . (4b) 

The coefficients a — e are complex but only the moduli \A\ and \B\ enter into the square 
brackets. We have determined all coefficients by comparing with bifurcation diagrams ob- 
tained for the flow at mid-gap position. 

We are seeking solutions of (jH) of the form A(t) = \A\e~ tQAt , B(t) = \B\e~ lUBt with 
constant moduli and frequencies. Here Qa,b = ^a,b—^° are the deviations of the frequencies 
from the critical ones according to (0 EH) • 

In view of Fig.Elwe have neglected here for simplicity reasons quintic contributions to SPI 
amplitudes and frequencies by discarding terms of the form c?|y4| 4 A and d\B\ 4 B in (j^,) and 
(jib), respectively. Taking a = 1 + one then has \A\ 2 = /j,/(—b r ) and Qa = — — bi\A\ 2 
for, say, a L-SPI (for a R-SPI replace A by B). The indices r and i denote real and imaginary 
parts, respectively. As an aside we mention that a linear variation of Qa with \A\ 2 holds for 
SPI also significantly further away from onset than in Fig. 121 

Motivated by Fig. El we also discard quintic contributions to the modulus and to the 
frequency of the RIB state. This is enforced by making in the terms /|I?| 4 v4 and f\A\ 4 B 
that appear on general grounds in (JUO and (jib), respectively, the special choice / = — e. 
Then the RIB solution is characterized by \A\ 2 = \B\ 2 = fi/(—b r — c r ) and Qa — &>b — 

-CLif! - (6j + Ci)\A\ 2 . 

The CR-SPI solution with finite A ^ B exists only for nonzero coupling e. It is most 
conveniently found and expressed in terms of the combined order parameters 

S =(\A\ 2 + \B\ 2 )/2 , D = (\A\ 2 -\B\ 2 )/2 (5a) 

n s = (n A + n B )/2 , n D = (n A -n B )/2. (5b) 

Then the radial velocity field (JHJ) of the CR-SPI at mid-gap can be written into the form u = 

^ A \ e i(kz-n n t) + |£| e -i(fcz-f2 t)j e i{Mip-u°t-n s t) + c c It showg the aforementioned rotation 

frequency us = uj° + Qs and the axial propagation of the deformation that is governed by 
the frequency = Qd- 
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The CR-SPI solution of ® reads 

S=(c r - 6 r )/(2e r ) , D 2 = + (b r + c r )S]/(2e r ) (6a) 
-n s = aiH+ (bi + Ci )S -2 ei D 2 , -tt D = (k - a + 2 ei S)D . (6b) 

Note that S is constant in this CR-SPI. D grows via a pitchfork bifurcation at /i2 = — (b r + 
c r )S' out of the RIB (D = 0) and connects at = —2b r S with D = \A\ 2 /2 = S with the 
L-SPI and with D = -\B\ 2 /2 = -S with the R-SPI. 

In Figs. 12]- HI we compare bifurcation diagrams of SPI, RIB, and CR-SPI solutions re- 
sulting from the amplitude equations (thin lines) with those from the full numerical 
simulations (symbols). For the sake of presentational clarity the stability of the former is 
not indicated in these figures. But stability analyses of the three vortex states in question 
confirms that the CR-SPI solution © transfers stability between the RIB and the SPI states 
as indicated in Fig. 

We should like to mention that keeping the discarded quintic terms by setting d ^ and 
/ 7^ — e slightly improves the agreement between the solutions of the Navier-Stokes equations 
and the amplitude equation approximation without changing the bifurcation topology and 
the stability behavior of the latter. However, then the CR-SPI solution would still have 
S = const. To obtain deviations from this behavior as displayed, e.g., by the triangles of 
Fig. Ufa) one needs higher orders in eqs. (jlj). 

So, the model (|IJ) is a minimal one in the sense that the quintic terms that are necessary 
to ensure the existence of A ^ B solutions enter only into the coupling between R- and 
L-SPI modes but not into the SPI and RIB solutions. Note furthermore that the coupled 
amplitude equations (jlj apply equally well to other systems with two distinct A and B 
solutions — be they stationary or oscillatory — that are related to each other by analogous 
(symmetry) requirements with similar nonlinear couplings. Thus, bifurcation diagrams in 
which the pure A, the pure B, and the A = B solutions are connected with a stability 
transfer by a A ^ B solution branch could also arise in these systems with a form that 
would be (topologically) similar to that of Fig. EJ 
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FIG. 1: (Color online) Snapshots of the radial flow u in the <p — z-plane of the 'unrolled' mid- 
gap cylindrical surface. White identifies the nodes. Light (dark) grey [red (blue) in the elec- 
tronic version] marks radial out-flow (in-flow). Arrows indicates how the phase moves. For 
the CR-SPI in (c) [(d)] the undulated pattern moves upwards [downwards] (full arrows) while 
the undulation itself propagates slowly downwards [upwards] (dashed arrows). Parameters are 
R 2 = -605.2(o, 6), -604.725(c, d), -604.5(e) and Rt = 240. 
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FIG. 2: (Color online) Bifurcation diagrams of SPI (red circles), RIB (blue diamonds), and CR- 
SPI (purple triangles) obtained from numerical solutions of the full Navier-Stokes equations for 
R\ = 240 (a,b) and R\ = 200 (c,d) as functions of R2 and /i. Shown are the squared mode am- 
plitudes \A\ 2 , \B\ 2 and frequencies ua,wb of the dominant modes 1*2,1(2), 1*2,-1 (t) ©, respectively, 
in the expansion (^Q) of the radial velocity field u at mid-gap. Full (open) symbols denote stable 
(unstable) solutions 12]. Lines show amplitude-equation approximations. Here a CR-SPI with 
\A\ > \B\,ua < transfers stability from an L-SPI (A ^ 0,B = 0) to the RIB state (A = B). 
The analogous diagram with an R-SPI is obtained by interchanging A and B. 
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FIG. 3: (Color online) Bifurcation diagrams of SPI (red circles), RIB (blue diamonds), and CR-SPI 
(purple triangles) obtained from numerical solutions of the full Navier-Stokes equations for R\ = 
240 using the order parameters (J5J): Relation between \x and S (a) and D (b), respectively; CR-SPI 
frequencies Qs = (^a + ^_b)/2 = {uja+^b) — w° (c) and = (^A~ ^s)/2 = {oja —ujr)/2 = u>d 
(d), respectively, versus D. Full (open) symbols denote stable (unstable) solutions 12]. Lines show 
amplitude-equation approximations. Here a CR-SPI (\A\ > \B\,D > and Qa < Qb,&d < 0) 
transfers stability from an L-SPI (B = 0, i.e., D = S) to the RIB state (A = B, i.e., D = 0). 



February 2, 2008 



11 




FIG. 4: (Color online) Bifurcation diagrams for R\ = 200 as in Fig. |21 
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FIG. 5: (Color online) Bifurcation diagrams of SPI, RIB, and CR-SPI resulting from the amplitude- 
equation approximation, (a) 3D plot in \A\ 2 — \B\ 2 — ^-space; (b) projection onto the \A\ 2 — (-Bp- 
plane or the D — (S-plane, respectively. Full (dashed) lines denote stable (unstable) solutions. The 
CR-SPI connects with constant S the two SPI solution branches to the RIB thereby transferring 
stability from the former to the latter. 



